Superconducting vortices in Semilocal Models 
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Abstract 

It is shown that the SU(2) semilocal model - the Abelian Higgs model with two complex 
scalars - admits a new class of stationary, straight string solutions carrying a persistent 
current and having finite energy per unit length. In the plane orthogonal to their direction 
fvq \ they correspond to a nontrivial deformation of the embedded Abrikosov-Nielsen-Olesen 

£> ■ (ANO) vortices by the current flowing through them. The new solutions bifurcate with the 

ANO vortices in the limit of vanishing current. They can be either static or stationary. 
In the stationary case the relative phase of the two scalars rotates at constant velocity, 
giving rise to an electric field and angular momentum, while the energy remains finite. 
The new static vortex solutions have lower energy than the ANO vortices and could be 
of considerable importance in various physical systems (from condensed matter to cosmic 
■ strings). 



It is difficult to overemphasize the importance of "topological defects" (monopoles, vortex 
lines, domain walls, etc) resulting from the spontaneous breaking of global or local symmetries 
in modern physics, we refer to the some of the numerous reviews [1]. There is an interest- 
ing class of models in which both global and local symmetries are simultaneously broken, in 
the most economical way, with the minimal number of scalar fields. These models, dubbed 
semilocal models, have some remarkable features, see the detailed review [2]. Perhaps the most 
remarkable point is that semilocal models exhibit vortices, despite the first homotopy group of 
the vacuum manifold being trivial [3, 4]. The simplest example of such a model is when electro- 
magnetism is coupled to two charged complex scalars transforming as a doublet under a global 
SU(2) symmetry (extended Abelian Higgs model). This simple case is made quite interesting 
by noting that it represents the bosonic sector of the Electroweak theory in the limit that the 
Weinberg angle, 0\y = 7r /2 [5]. Furthermore such multicomponent Ginzburg-Landau models 
are also used in condensed matter theory, for example in the description of unconventional 
superconductors [6], in modeling neutral two-component plasmas [7]. The Abrikosov-Nielsen- 
Olesen (ANO) vortex solutions [8], characterized by the mass ratio, (3 = m s /m v (m s resp. m v 
denoting the mass of the scalar resp. of the vector fields), are embedded into this theory. The 
parameter f3 also distinguishes superconductors of type I (f3 < 1) from type II (/? > 1). The 
case (3 = 1 is quite special for here, instead of there being a unique vortex solution (the ANO 
vortex), there is a family of them labelled by a complex number of the same energy, satisfying 
the Bogomol'nyi equations [9]. The stability of the embedded ANO vortices in Abelian models 
with an extended Higgs sector was examined in Ref. [4], and it was shown there that they are 
stable only if (3 < 1. The counterparts of these solutions are the Z-strings in the Standard 
model which have been shown to be stable for sin 2 Q-yj > 0.9 [2]. 

In this Letter we present a new class of straight flux-tube solutions (strings) of the SU(2) 
symmetric semilocal model which are translationally symmetric along the z = £3 direction. 
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The new solutions form a one parameter family for any fixed value of (3 > 1. The main feature 
distinguishing them from the previously known string or equivalently vortex solutions in the 
(xi,X2) plane, is due to the presence of a relative phase between the two scalars, depending 
linearly on z. Because of this relative phase (twist) there is also a persistent current flowing 
along the z-direction, so in fact these new objects are not genuinely planar vortices. The 
relative phase may also have a linear dependence on the time coordinate, xo, in which case the 
vortices are stationary and possess a nonzero momentum, angular momentum, as well as an 
electric field. Due to the strong (exponential) localization of the fields the vortices have finite 
energy per unit length. The energy of the current carrying solutions is lower than that of the 
corresponding ANO vortex, therefore they should have better stability properties. When the 
current vanishes the solution bifurcates with the embedded ANO vortex. In fact the well known 
"magnetic spreading" instability of the embedded type II ((3 > 1) ANO vortex in the semilocal 
model can be understood as a signal of bifurcation with a new solution branch, which carry a 
nonvanishing current. The Z-string counterparts of the new solutions should also be present in 
the Electroweak model and they may have better stability properties because of the localizing 
effect of the current. Our solutions can also be interpreted as superconducting cosmic strings 
of the bosonic type, but without an additional U(l) gauge field as in the prototype U(l)xU(l) 
model [10]. There are also some important differences between the superconducting strings in 
Witten's model as compared to the semilocal ones. For example the current carrying semilocal 
strings are strongly localized. 

The (suitably rescaled) Lagrangian of the SU(2) semilocal theory can be written as 

c = h {-\ F ^ F " V + lD ^ 2 - T m2 - 1)2 } (1) 

where <3? = (<pi, ^2), = d^ — iA^, |$| 2 = (ftifii + 4>2<p2, and g is the coupling constant, which 
will be set to 1 for simplicity. The breaking of the U(l) gauge symmetry leads to the physical 
spectrum consisting of a vector particle with mass m v = gr/, two Nambu-Goldstone bosons, 
and a Higgs scalar with mass m s = f3r), rj is the vacuum expectation value of the scalar fields. 
The most general translationally symmetric (in the X3 = z direction) and stationary Ansatz 
can be written as 

A,j, = (A (x 1 ,x 2 ) , A i (x 1 ,x 2 )) i = 1,2,3 
0i = /i(xi,x 2 ) 2 = / 2 (x 1 ,x 2 )e^" t +^), (2) 

where (coo, LU3) = w a , are real parameters. The most important feature of the Ansatz (2) is that 
the change in the phase of the scalar fields under an infinitesimal translation is compensated 
by a nontrivial gauge transformation. In other words a space-time translation moves the 
fields along gauge orbits. Also the Ansatz (2) reduces the global SU(2) symmetry to U(l). 
The Noether current associated to the remaining global U(l) symmetry is then simply = 
i(0 2 D At 2 — (t>2D^(t> 2 )/2. Therefore the solutions described by the Ansatz (2) have a conserved 
Noether charge, 1q, as well as a current, J3, flowing along the string where 

l a = ±r] 2 J d 2 x{u a - A a )4> 2 4>2 , a = (0,3). (3) 

The equations of motion imply the following "Gauss-law" constraints on A a = (Aq, A3) in the 
(xi, X2) plane: 

AA a - A a \$\ 2 + uj a j)2<i>2 = 0. (4) 

Now exploiting the "constraint" equations (4) one can establish the following lower bound on 
the total energy per unit length: 

E > ir\n\r] 2 + \uJ a 2 a \ , (5) 
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where n is an integer determining the magnetic flux. Furthermore from Eqs. (4) it also follows 
that for any regular solution Aq = ujqA 3 /uj 3 [uj 3 7^ 0), i.e. Aq is determined by A 3 . This is 
related to the fact that one can still perform Lorentz transformations (boosts) in the (t, z) 
plane, which mix the Aq and A3 components of the vector potential, and also transform the 
parameters (ujq ,uj 3 ) of the Ansatz (2). In fact it is only their Lorentz invariant combination, 
uj 2 = uj 2 — uj 2 , that appears in the equations of motion. Therefore if to 2 > 0, by a Lorentz boost 
one can always achieve ujq = , Aq = 0, i.e. it is sufficient to consider the static case only. It 
should be pointed out here that the general Ansatz (2) has been also considered by Abraham 
[11] in the same model. He considered, however, exclusively the Bogomol'nyi equations existing 
only in the very special case (5=1. As found in Ref. [11] this implies further uj 3 = ujq i.e. uj 2 = 0, 
which is referred to as the chiral case, when due to the lack of localization the energy of the 
solutions exhibited in Ref. [11] (with winding number n = 1) diverges. 

To simplify further the problem of finding current carrying string solutions of the semilocal 
theory Eq. (1), we also impose rotational symmetry in the (x\,X2) plane, leading finally to the 
following Ansatz in polar coordinates: 



Aq = a (r) , A r = , A^ = na(r) , A 3 = uja 3 (r) , 
fa = / 1 (r)e i "^ fa = / 2 ( r y™/v(^>i+^32) j 



(6) 



where m = 0, ...n— 1. In this Letter we shall concentrate on the simplest m = case, a 
more exhaustive investigation will be published [12]. The static (i.e. ao = 0) field equations for 
m = become: 
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ra' 3 )' = a 3 \f\ 2 -f 2 , with >=d/dr, 



rC-)' = a\f\ 2 -f 2 



(r/D' = /i 
(rf 2 )' = f2 



2 (1 a ) 1 '• 2 ' , r 1 /• |2 , 
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n 
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(7a) 
(7b) 
(7c) 

(7d) 



Eqs. (7) admit a 4-parameter family of local solutions regular at the origin where a = a^r 2 + 



0{r%a 3 = af+0{r 2 ),h 



f^r n + 0(r 2+n ), h = f2 } +0(r 2 ) where , af, tf\ / 2 (0) are 
free parameters. For r — ► 00 the regular solutions decay exponentially: a = 1 + O (exp [— br]), 
a 3 = O (exp [— br]), /1 = 1 + O (exp [— cr]), f 2 = O (exp [— ur] ) , where the constants, b, c, 
determining the exponential decay are given as b = min(l,2u;), c = min(2, f3, 2u), and they 
also comprise a 4-parameter family of local solutions. From this asymptotic behaviour it is 
clear that uj > is an essential parameter to ensure exponential localization. The reduced 
energy functional takes the form: 



E 
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rdr { n z ^ + u?a 3 2 + \j ri \ + n 



p/2i 



(1 - a)' 



ft + 3/1 



^ 2 [/i 2 ^ + /l(l-a3) 2 ]+^(l-|/ 2 |) 2 }, 



+ 



(8) 



from which it is apparent that the exponentially localized solutions have finite energy. In view 
of the mathematical complexity of Eqs. (7) we have resorted to numerical techniques to obtain 
regular solutions. Adapting a multishooting procedure [13], we have succeeded in obtaining 
numerically families of regular solutions of Eqs. (7), representing a new class of vortices; see 
Fig. 1 for a sample vortex profile. For a given value of (3 > 1 our solution class is parameterized 
by uj, which varies in a finite interval < uj < uJb(P)- As the parameter uj tends to its maximal 
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Figure 1: A superconducting vortex solution for (3 = 2,oj = 0.1. 



value, u> — > u\>{(3), the functions 03 — ► 0, fa — ► 0, while a — > a ANO , /1 — ► /ano> i- e - the new 
class bifurcates with the corresponding ANO vortex. This bifurcation can be understood by 
linearizing Eqs. (7) around the ANO solution with 03 « 1, /2 « 1. To first order 03 = and 
the linearized fa amplitude satisfies 



--(rfa)' + 
r 



IM) _ P /-I _ r 2 

2 o V 1 JanoJ 



fa = -ul((3)fa. (9) 



Eq. (9) corresponds precisely to the linearized stability problem of the embedded ANO vortex 
in the semilocal theory studied by Hindmarsh [4]. In agreement with his results 1 , for any 
(3 > 1 we find a single normalizable solution, fa, with negative eigenvalue, — In Ref. 
[4] this bound state is a sign of instability of the semilocal ANO vortex. In our case the very 
same bound state indicates the bifurcation of a new branch with that of the ANO solution at 
co = Ub(fl). Quite importantly the energy of the new solutions (with n = 1) is always smaller 
than that of the ANO vortex. Therefore we expect the new solutions of winding number n = 1 
with a nonzero current to be stable. 

It is instructive to discuss also the (3 = 00 limiting theory, when the scalar fields are 
constrained |/i| 2 + \ fa\ 2 = 1, (a gauged CP 1 -model) whose field equations can be written as: 

i(ra' 3 )' = o 3 -sin 2 fl, r(^)' = a - cos 2 9, (10a) 



uj 2 (1 - 2a 3 ) - n 2 



sin(20) , (10b) 



where we have introduced fa = cos 9, fa = sin 9. The profile functions of the current carrying 
vortex solution of Eqs. (10) are qualitatively similar to those of Eqs. (7) (comp. Fig. 1). The 
energy of the current carrying solutions is finite for any value of uj, whereas the energy of the 
[3 = 00 ANO vortex diverges. This difference is also reflected by the fact that ^b(°°) = 00, i.e. 
the parameter to varies in the interval (0, 00). 

The phase space of the new solution class is illustrated on Fig. 2, choosing the total 
energy, E, the magnitude of the condensate at the origin, and uj as parameters. On Fig. 
2 a family of solutions for a fixed value of j3 corresponds to a curve in the (f^ , oj) plane, and 

1 Note that our (3 2 corresponds to f3 of Ref. [4]. 
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Figure 2: Phase space of the n = 1 superconducting vortices. 



we have indicated the corresponding ANO vortex, located at the endpoint of the curve in the 
= plane where the new family bifurcates with it. The family of solutions for (3 = oo 
is represented by the curve in the = 1 plane. As the value of oj decreases from uJb(P) 
towards oj = 0, the energy of all superconducting solutions decreases towards the Bogomol'nyi 
bound, E = nr) 2 , without ever attaining it, irrespectively of the value of (3. An important point 
is that solutions with nonzero values of oj do not have a smooth oj — > limit, i.e. the point 
(fi^ = 1)^ = 0) does not belong to the phase space as it is indicated on Fig. 2 by a circle. In 
the (5 = oo theory one can obtain the following virial relation 

j\dr |n 2 ^-i,j 2 [a'3 + (l-2a3)sin 2 0]| =0, (11) 

which immediately implies that for oj = the (3 = oo theory does not admit nontrivial finite 
energy solutions. For finite values of /3, a suitable change of scale in Eqs. (7) shows that the 
effective value of the scalar self coupling is given by j3 c g = (3/uj. It follows that the limit oj — ► 
also implies f3 c g — > oo, enforcing |/i| 2 + I/2I 2 —> 1. This explains the universal behaviour of the 
solutions in the oj — > limit. A more detailed analysis shows [12] that all solutions converge 
pointwise to a limiting configuration corresponding the CP 1 lump analyzed in Ref. [4]. 

In the limit oj — > the value of the current, X3, is getting large and it even seems to 
diverge while 0JI3 — > 0. This behaviour of the current is in sharp contrast with what is found in 
the U(l) xU(l) model, where the current has a maximal value above which the superconducting 
string "goes normal" [14]. The z-component of the current, X3 = Iz/ftirr] 2 ), is depicted on Fig. 
3 as a function of the parameter oj. 

Since our vortex solutions exist only in type II superconductors (j3 > 1) it is natural 
to ask as to their interaction energy (which is well known to be repulsive for the type II 
ANO vortices). Since we have superimposed vortex solutions (for n > 1) we can ask if it 
is energetically favourable for them to break up into widely separated constituents. To this 
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Figure 3: The current, T 3 as a function of lo for n = 1, f3 = \/2, 2, 3, oo. (The curves terminate 
at the smallest value of to we could compute) 
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Figure 4: Energy /winding, E/^nrj 2 ), of the n = 2,3 vortices for (3 = y/2 in function of the 
current, Z3. 

end we depict E{n,T%) / {imrf) , the energy per winding number of the superimposed n = 2 and 
n = 3 solutions as a function of the current flowing through them on Fig. 4. The graphs on Fig. 
4 indicate that the interaction is repulsive up to a certain value of T3, however, at least for the 
n = 2 vortex it becomes attractive when X3 > 0.71. This suggest that the n = 2 superimposed 
vortex for large enough currents becomes stable with respect to breakup into n = 1 vortices 
unlike the type II ANO vortices. If the interaction for separated vortices has also an attractive 
phase this could have important physical consequences (e.g. the attraction may significatively 
change the intercommutation of colliding strings). Let us mention here that in a somewhat 
different two-component Ginzburg-Landau model similar non-monotonic behaviour of E/n has 
been reported in [15]. 

Stationary solutions (with ujq / 0) are obtained by a Lorentz boost of the static ones along 
the z-axis. The stationary vortices obtained this way have momentum, P = loqI^/lo^, angular 
momentum, J = —ujQnls/uj^, and they are surrounded by a radial electric field E r = ujQa' 3 . 
It is interesting that in our case the electric field is screened, and the total energy of the 
superconducting vortices stays finite. 
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In conclusion such current carrying strings provide a new interesting class of defects 
which may be realized in realistic physical systems described by semilocal models and they 
may equally find their applications in a cosmological context as cosmic strings and may even 
be present in the Standard model of Electroweak interactions. Finally we thank G. Volovik for 
pointing out that twisted vortices have also been investigated and even experimentally observed 
in superfmid 3 He [16]. 
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